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Abstract 

The molecular diversity and community structure of bacteria from the gut of the termite Reticulitermes speratus were analyzed by the 
sequencing of near-full-length 16S rRNA genes, amplified by polymerase chain reaction. The results of the analysis of 1344 clones 
indicated a predominance of spirochetes in the gut. Spirochetal clones accounted for approximately half of the analyzed clones. The clones 
related to Bcicteroides, Clostridia, and the candidate division Termite Group I each accounted for approximately 5-15% of the analyzed 
clones. The rest were comprised of Proteobacteria, Actinobacteria, Mycoplasma and others. Using the criterion of 97% sequence identity, 
the clones were sorted into 268 phylotypes, including 100 clostridial, 61 spirochetal and 31 Bcicteroides - related phylotypes. More than 90% 
of the phylotypes were found for the first time, and some constituted monophyletic clusters with sequences recovered from the gut of 
other termite species. 

© 2003 Federation of European Microbiological Societies. Published by Elsevier Science B.V. All rights reserved. 
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1. Introduction 

Termites harbor an abundance and diversity of gut mi- 
crobes, comprised of bacteria, protists, fungi and archaea. 
The termite gut microbe ecosystem is one of the most 
fascinating examples of symbiosis between an animal 
and microbes, and among a diversity of microbes. The 
characteristics of these gut microbes have been extensively 
studied, and it is believed that the microbes provide car- 
bon, nitrogen and energy nutrition to their host termites, 
to the extent that termites can no longer live without them 
[1-3]. However, difficulties in the cultivation of these gut 
microbes and the complexity of their community structure 
have been obstacles to a comprehensive understanding of 
the symbiotic system. To overcome these problems, efforts 
to cultivate and characterize as yet unculturable microbes 
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have continued [4,5] and culture-independent approaches 
have been applied [6-11]. In our previous work, partial 
16S rRNA genes were amplified by polymerase chain re- 
action (PCR) from the gut of the termite Reticulitermes 
speratus to estimate the diversity of the gut bacterial mi- 
crobiota [6], We found a predominance of five bacterial 
groups: the Cytophaga-Flexibacter-Bacteroides (CFB) 

group; the low-G+C Gram-positive bacteria; Proteobac- 
teria; Spirochaeta; and the newly discovered ‘Termite 
Group I’ (TG I). However, we analyzed only 55 clones 
and sequence lengths of approximately 300 bp in that 
study, insufficient for the reliable elucidation of bacterial 
microbiota and their phylogenetic affiliation. While the 
investigation of microbiota using PCR amplification and 
the sequencing of small subunit rRNA genes contains 
some problems, such as PCR errors, PCR bias and differ- 
ent copy numbers of the gene, it is widely accepted as a 
powerful and essential tool for the investigation of as yet 
unculturable microbial communities. In studies of termite 
gut microbe symbiosis, thorough molecularly based inves- 
tigations of microbiota are needed not only for studies of 
their phylogenetic diversity and evolution, but also for 
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evaluating the ecological and physiological impact of var- 
ious gut microbes on a termite’s metabolism. 

In the present paper, we report a detailed analysis of 
bacterial diversity and community structure in the gut of 
the termite R. speratus. Data are based on sequences of 
near-full-length 16S rRNA genes, amplified by PCR. To 
reduce the possibility of overlooking any major bacterial 
populations in the analysis, several sets of Bacteria-univer- 
sal primers and annealing temperatures were used. 

2. Materials and methods 

2.1. Sample collection and extraction of DN A 

Wood-feeding lower termites, R. speratus (family Rhi- 
notermitidae), were collected at Ogose, Saitama, Japan. 
Within 3 h of collection, the termites were washed in steri- 
lized water and their guts were drawn out using sterilized 
forceps. The isolated guts of 30 individuals were sub- 
merged in solution U [11] on ice and gently crushed using 
a pestle. After centrifugation at 2000 Xg for 5 min, the 
precipitate was treated with 0.2 mg ml -1 Zymolyase 
100-T (Seikagaku) and 2 mg ml -1 lysozyme (Nacalai) at 
30°C for 30 min. The reaction mixture was dropped into 
liquid nitrogen and crushed into a fine powder, from 
which DNA was extracted using the Qiagen RNA/DNA 
kit. 

2.2. PCR amplification and preparation of clone libraries 

16S rRNA genes were amplified from the gut sample by 
PCR, using Bacteria-universal primers (Table 1). PCR was 
performed using a PTC-200 Thermal Cycler (MJ Re- 
search), EX-Taq polymerase (Takara), and the following 
program: 2 min of initial denaturation at 95°C followed 
by 12 or 18 cycles of denaturation (30 s at 95°C), anneal- 
ing (60 s at 45, 50 or 55°C), and extension (4 min at 72°C), 
with a final extension at 72°C for 10 min. The concentra- 
tions of the template and primers were 0.1 ng pi -1 and 
1 pM, respectively. Fourteen sets of PCRs, using different 
primers and the above annealing temperatures and cycle 
numbers, were performed to prepare independent clone 


Table 1 

Bacteria-universal primers for PCR and sequencing 


Primer 

Sequence 

Reference 

27F 

5'-agagtttgatcmtggctcag-3' 

[38] 

39F 

5'-tggctcagrwygaacgctrg-3' 

This study 

41F 

5'-gctcagattgaacgctggcg-3' 

[39] 

63F 

5'-caggcctaacacatgcaagtc-3' 

[40] 

64F 

5'-bgyctwanrcatgcaagtyg-3' 

This study 

533F 

5 '-tgccagcagccgcggtaa-3 ' 

[41]“ 

908R 

5'-cgtcaattcmtttgagtt-3' 

[41] 

1389R 

5'-acgggcggtgtgtacaag-3' 

[40,42] 

1492R 

5'-ggttaccttgttacgactt-3' 

[38] 


a Slightly modified from the original primer. 


libraries (Fig. 1). In each PCR set, the products were pu- 
rified from 200 pi of reaction mixture, using the QIAquick 
PCR Purification Kit (Qiagen), followed by ethanol pre- 
cipitation. The products were cloned into pCR-4 TOPO 
TA cloning vectors (Invitrogen), and clone libraries were 
established on LB agar plates. 

2.3. Sequencing and phylogenetic affiliation of clones 

Ninety-six clones were randomly chosen from each 
clone library for analysis. Approximately 550 bp of the 
16S rDNA clones, corresponding to the region 563-1114 
in Escherichia coli (J01695), were sequenced with primer 
533F (Table 1), using the Big-Dye Terminator Cycle Se- 
quencing Kit (Perkin-Elmer) and ABI 3700 or 377 genetic 
analyzers. 

The clones sharing more than 98% sequence identity 
were grouped into the same phylotype, and the full length 
of representative clones of all the phylotypes was se- 
quenced with primers T7, T3 and 908R (Table 1), using 
the long read protocol distributed by PE Biosystems. They 
were further sorted into phylotypes using the criterion of 
97% sequence identity. All the sequences were subjected to 
the Check-Chimera program on the Ribosomal Database 
Project (RDP) web site [12] for the elimination of chimeric 
sequences, and analyzed by BLAST 2.0 [13] and FASTA 
3.0 [14] to find the closest phylogenetic neighbors. The 
sequences of all the phylotypes were aligned with those 
which marked the highest point in the FASTA analysis 
and with other representative sequences from on-line data- 
bases, using CLUSTAL X1.8 [15]. The region correspond- 
ing to positions 64-1388 in E. coli (JO 1695) was used for 
the alignment, and all the sites containing gaps were ex- 
cluded from the phylogenetic analysis to decrease ambigu- 
ity. As for shorter sequences, only those which were 
thought to be essential were included, and the sites con- 
taining missing data were excluded pairwise in the phylo- 
genetic analysis. A neighbor-joining tree of the aligned 
sequences was constructed using MEGA V2.1 [16], based 
on the distance matrix calculated according to the Jukes- 
Cantor model. The tree was tested statistically using a 
bootstrap resampling of 1000 times. GENETYX V10.1.6 
(Software Development) was used to calculate sequence 
similarities. DNASIS V3.7 (Hitachi) was used for the gen- 
eral manipulation of the sequence data. The sequence data 
will appear in the DDBJ, EMBL and GenBank nucleotide 
sequence databases under accession numbers AB088857- 
AB089126. 

2.4. Statistical analysis 

The difference between clone libraries was examined 
pairwise by the statistical analysis described by Singleton 
et al. [17], From each clone library, 96 sequences of ap- 
proximately 550 bp, corresponding to the region 563-1114 
in E. coli (J01695), were used for inter-clone comparisons. 
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Fig. 1. Relative abundance of 16S rDNA clones in the 14 clone libraries A-Q. The primer sets, annealing temperatures and number of cycles for PCR 
are attached to the library names. Ninety-six clones were analyzed for each clone library. The clones were affiliated with the taxonomic groups by phy- 
logenetic analysis. Rare taxa were grouped as ‘others’. 


Distance matrices were calculated using MEGA V2.1, ac- 
cording to the Jukes-Cantor model and based on the se- 
quences aligned by CLUSTAL XI. 8, and were analyzed 
using the program LIBSHUFF [17], distributed freely at 
http ://www. arches. uga.edu/ ~ whitman/libshuff.html. The 
coverage [18] was calculated by the formula [1— (nIN)], 
where n is the number of phylotypes represented by only 
one clone and N is the total number of clones. 


3. Results and discussion 

3.1. Estimation of bacterial diversity and community 
structure in the termite gut 

From the total of 1344 (96 X 14) clones, 268 phylotypes 
were recognized. The phylotypes were classified phyloge- 
netically into 1 1 bacterial divisions, an unidentified cluster 
and an unidentified sequence (Table 2). The relative abun- 
dance of clones affiliated with the respective taxonomic 
groups (Table 2) is shown in Fig. 1. No statistically sig- 
nificant difference between the clone libraries was detected 
using the program FIBSFIUFF [17], except in the compar- 
ison between a library prepared from PCR using primer 
63F and a library from PCR using the other primers 
(■ P <0.01). The difference was obviously due to the lack 
of spirochetal clones in libraries K, M and N (Fig. 1), 
presumably caused by preferential annealing due to six 
or seven mismatched base pairs between 63F and the cor- 
responding region of the spirochetal sequences. These 
three clone libraries were not used in the estimation of 
community structure and diversity. 


In all the clone libraries, except for K, M and N, spi- 
rochetal clones were the most dominant, accounting for 
42-63% of the sequenced clones (Fig. 1). This suggests a 
predominance of spirochetes in the gut of R. speratus, 
which is comparable to the result of the microscopic ob- 
servation of the gut contents of the higher termite Nasu- 
titermes lujae [19]. The Bacteroides-related, clostridial and 
TG I clones formed the second most dominant groups, 
accounting for 3-19%, 6-16%, and 4-11% of the se- 
quenced clones, respectively. 

The community structure at the phylotype level was out- 
lined in the species abundance curve depicted in Fig. 2. 



Fig. 2. Species abundance curve for the 16S rDNA clone libraries A-J 
and P-Q. The data from the clone libraries K, M and N are excluded 
because of their lack of the most dominant, spirochetal clones. The total 
numbers of clones and phylotypes were 1056 and 218, respectively. The 
16S rDNA phylotype Rs-D73 is indicated as the most abundant. 
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Table 2 


Classification and number of 16S rDNA phylotypes found in this study 


Taxon 

RDP code 

Number of phylotypes 1 * 

Accession number 

Spirochaeta 




Treponema 

[2.27.3.2] 

57 

AB088857-916 

Spirochaeta 

[2.27.3.1,3] 

3 


Undescribed cluster 

none b 

1 

AB089126 

CFB group 




Bacteroides and relatives 

[2.15.1] 

31 

AB088917-47 

Low-G+C Gram-positives 




Clostridium! Eubacterium 

[2.30.4, 5, 9] 

100 

AB088948-9047 

Mycoplasma 

[2.30.8.4] 

7 

AB089053-60 

Lactococcus, Enterococcus 

[2.30.7.20, 21] 

3 

AB089061-3 

Anaerobaculum 

[2.11] 

3 

AB089064— 6 

Termite Group I 

[2.19] 

5 

AB089048-52 

Actinobacteria 

[2.30.1] 

14 

AB089069-82 

Proteobacteria 




a 

[2.28.1] 

16 

AB089083-99 

p 

[2.28.2] 

4 

AB089100-3 

5 

[2.28.4.1. 10] 

6 

AB089104— 10 

8 

[2.28.5] 

3 

AB0891 1 1-3 

Y 

[2.28.3] 

2 

AB0891 15-6 

Planctomycetes 

[2.20] 

4 

AB0891 17-20 

V errucomicrobia 

[2.10] 

2 

AB089121-22 

Cyanobacteria 

[2.21] 

1 

AB089123 

Acidobacteria 

[2.25.3.6] 

1 

AB089124 

OPll 

[2.28.4.2] 

2 

AB0891 14, -125 

Unidentified cluster 

none 

2 

AB089067-8 

Unidentified bacterium 

none 

1 

AB089107 

Total number of phylotypes 


268 



a Phylotypes were defined using the criterion of 97% sequence identity in the region corresponding to positions 64—1388 in E. coli (J01695). 
b The reference sequences from the databases in this cluster are classified into different RDP taxonomic groups as follows. SJA-88: Nitrospina subdivi- 
sion [2.13.6]; vadinBA43: Heliobacterium group [2.30.3.2] in the low-G+C Gram-positive bacteria; EKHO-6; Thermoanaerobacter and relatives [2.30.2] 
in the low-G+C Gram-positive bacteria. Also see Fig. 3. 


Whereas a species abundance of prokaryotes is theoreti- 
cally deduced to approximate a log-normal distribution 
[20], namely, few numerically dominant and few very 
rare species but many species of intermediate abundance, 
as demonstrated in the gut bacteria of pigs [21], it was 
quite different in the present study. We found few domi- 
nant phylotypes and many rare phylotypes (Fig. 2). The 
number of sequenced clones in the study of gut bacteria 
from pigs [21] was up to 4270, much more than the 1056 
(without the clones from libraries K, M and N) in our 
study and the coverage, 0.978 [21], was also much higher 
than 0.905 in the present study, i.e., 2.2% of the addition- 
ally sequenced clones from pig intestine would represent 
new phylotypes while 9.5% of those in the present study 
would. These differences imply that the number of se- 
quenced clones in the present study is insufficient to reflect 
the diversity of all bacteria in the termite gut. However, 
Curtis et al. [20] recently proposed a mathematical equa- 
tion to estimate prokaryotic diversity in an environment 
from a practically acquirable data set of 16S rDNA clone 
libraries. To solve the formula [20], we need the values Nj, 
A max and A m ; n , where At is the total number of individ- 
uals in the community, A max is the number of individuals 
in the most abundant species and A m ; n is the number of 


individuals in the least abundant species. Aj can be esti- 
mated from a previous study [22] using Reticulitenn.es fla- 
vipes, a congeneric of R. speratus with similar body size. In 
R. flavipes, the total count of prokaryotic cells was around 
10 10 cells mf 1 gut fluid, corresponding to approximately 
3X10 6 cells gut -1 [22]. The population of archaeal cells 
can be ignored because they were estimated to occupy less 
than 10% of the prokaryotic community [10]. A max can be 
estimated from At and the frequency of the most abun- 
dant phylotype, Rs-D73, which comprised approximately 
10% of the sequenced clones from all the clone libraries, 
except for K, M and N. Assuming A m ;„ is 1, as suggested 
by the authors [20], the equation was solved for the num- 
ber of phylotypes to be about 6000 phylotypes ml -1 gut 
fluid and 700 phylotypes gut -1 . However, we should be 
careful in handling this type of approximation. Besides the 
intrinsic problems of this estimation, especially when a 16S 
rDNA clone library is used as the data [23], using the 
frequency of Rs-D73 as the most abundant phylotype 
could be problematic. Since this phylotype was identical 
to one of the dominant bacterial phylotypes amplified by 
PCR from the gut protist Dinenympha porteri [11], it is 
possible that it is one of the protist ectosymbionts 
and that its habitat is restricted to the cell surface of 
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D. porteri. In this case, the abundance of Rs-D73 might 
not be appropriate for estimating diversity in the whole 
gut. To more accurately estimate the microbial diversity in 
termite guts, further ecological studies, including the local- 



ization and characterization of dominant bacteria, are 
needed. 

3.2. Phylogenetic affiliation of 16S rDNA phylotypes 
3.2.1. Spirochaeta 

As described above, spirochetes comprised approxi- 
mately half of the sequenced clones. Sixty-one phylotypes 
were affiliated with this group and 51 of them were found 
for the first time. Most of the phylotypes (57 of 61) were 
affiliated with the genus Treponema (Table 2 and Fig. 3). 
They were further classified into two distinct clusters, ‘ter- 
mite Treponema clusters I and IF, as designated in our 
previous studies [11,24]. Fifty phylotypes were affiliated 
with cluster I and seven with cluster II. The sequence 
similarities between the 50 phylotypes in cluster I were 
more than 87% and those of the seven phylotypes in clus- 
ter II were more than 88%. The most abundant phylotype 
in all of the sequenced clones, Rs-D73, was included in 
cluster I and was almost identical ( > 99%) to RsDiSp8 
(AB032006), which was recovered from the protist D. por- 
teri in the gut of R. speratus [11], The second most abun- 
dant phylotype in all of the sequenced clones, Rs-G65, was 
affiliated with cluster II and almost identical ( > 99%) to 
RsDiSp9 (AB032002), also recovered from D. porteri 
[11], The phylotype Rs-EOl, sharing 98.7% sequence iden- 
tity with the D. porteri- derived sequence RsDiSpl (AB- 
031997), was also abundantly found in the sequenced 
clones (Figs. 3 and 4). These phylotypes implied an abun- 
dance of spirochetal ectosymbionts with the protists in the 
gut of R. speratus. In addition to these three phylotypes, 
several abundant phylotypes were included in the Trepo- 
nema clusters as shown in Fig. 4. 

Although the physiological characteristics of bacteria 
cannot be deduced from their phylogenetic position alone, 
the predominance of treponemes in the gut of R. speratus, 
as suggested in this study, implies interesting possibilities. 
The Treponema spp. in cluster I, ZAS-1 and ZAS-2, iso- 
lated from the gut of the termite Zootermopsis angusticollis 
were characterized as H 2 /CO 2 acetogens [4], and ZAS-9 
was characterized as a N 2 fixer [5], These microbial activ- 
ities are thought to be essential for the nutrition of wood- 
feeding termites [2], Therefore, it is worth attempting to 


Fig. 3. Phylogenetic tree showing the positions of 16S rDNA phylo- 
types affiliated with Spirochaeta, recovered from the gut of the termite 
R. speratus. The total number of clones belonging to each phylotype is 
indicated in parentheses. Only bootstrap confidence values greater than 
50 are indicated. The scale bar represents a 3% estimated sequence di- 
vergence. The ‘undescribed cluster’ includes the short sequences MERTZ 
2CM-174 and EKHO-6 (corresponding to positions 518-1388 in E. coli 
JO 1695). EKHO-6 was found in the databases as a member of the can- 
didate division Sediment-4. The origins of the reference sequences from 
termite guts are as follows. za29: Zootermopsis angusticollis ; RsDiSp se- 
ries: R. speratus [11]; HsDiSp320 and HsPySpl5: Hodotermopsis sjoes- 
tedti [11]; RFS series: R. flavipes [9]; ZAS series: Z. angusticollis [5,9]; 
nc26: Neotermes castaneus; BCf series: Coptotermes formosanus. 
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Fig. 4. Abundance of 16S rDNA phylotypes frequently found from all 
the sequenced clones. Since the clone libraries K, M and N lacked spi- 
rochetal clones, their data are distinguished from the others. Only phy- 
lotypes represented by more than 10 clones are shown. 

characterize abundant phylotypes such as Rs-D73, Rs- 
G65 and Rs-A48 for these activities using culture-depen- 
dent and/or -independent methods. 

Only a few spirochetal phylotypes were affiliated outside 
Treponema. Three phylotypes were affiliated with the ge- 
nus Spirochaeta and two of the three constituted a mono- 
phyletic cluster with za29 (AJ419823), the sequence ampli- 
fied from the gut of Z. angusticollis (Fig. 3). The phylotype 
Rs-H88 was clustered with four unidentified environmen- 
tal samples and its taxonomic position is uncertain (Table 
2 and Fig. 3). 

3.2.2. Bacteroides and relatives 

The Bacteroides- related clones were sorted to 31 phylo- 
types; 28 of them were found for the first time. The three 
phylotypes Rs-N74, Rs-B92 and Rs-M21 were almost 
identical (>99%) to the sequences UN71 (D63629), 
UN78 (D63632) and UN108 (AB055736), respectively, 
which were also recovered from the gut of R. speratus 
[6], Since these ‘UN-series’ sequences were as short as 
about 850 bp, we excluded them from the phylogenetic 
analysis. One of the three, Rs-N74, was the most abun- 
dant Bacteroides-related phylotype from the sequenced 
clones (Figs. 4 and 5). It constituted a monophyletic clus- 
ter exclusively comprised of termite gut-derived clones 
(corresponding to cluster V in our previous study [25]), 
although the bootstrap confidence values for the roots 
were not significantly high (Fig. 5 and [25]). The sequences 
listed on the top of the BLAST results for 19 Bacteroides- 
related phylotypes, including the above three, were from 
the guts of various termites. Since all the related termite- 
derived sequences were as short as about 850 bp except for 
the BCf series from the termite Coptotermes formosanus 
(only found in databases as the clone series ‘BCf), they 


were not included in the present analysis. Instead, a de- 
tailed phylogenetic analysis of these shorter Bacteroides- 
related sequences, including the above UN series, was 
done and reported in our previous paper [25]. The se- 
quence similarities between the 3 1 Bacteroides-related phy- 
lotypes were 78-95%. 

3.2.3. The low-G+C Gram-positive bacteria 
The clones affiliated with the low-G+C Gram-positive 
bacteria were very diverse and sorted to 113 phylotypes. 
All of them were found for the first time. Most of them, 
100 of 113, were affiliated with Clostridia (Table 2), con- 
sisting of 44 phylotypes related to Clostridium leptum and 
C. thermocellum (Fig. 6), 23 phylotypes related to bovine 
ruminal clones, swine intestinal clones and termite gut 
clones in unidentified clusters (Fig. 7), 23 phylotypes re- 
lated to C. propionicum and C. xylanolyticum (Fig. 8A), 



Fig. 5. Phylogenetic tree showing the positions of 16S rDNA phylotypes 
related to Bacteroides in the CFB group, recovered from the gut of the 
termite R. speratus. See the legend of Fig. 3 for explanation. 
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and 10 phylotypes related to C. purinolyticum (Fig. 8B). 
All of these phylotypes were found in small numbers (one 
to nine) from the sequenced clones. The BLAST analysis 
of the 100 clostridial phylotypes revealed that the database 
sequences showing the highest homology to 28 phylotypes 



ue 


Fig. 6. Phylogenetic tree showing the positions of 16S rDNA phylotypes 
related to C. leptum and C. thermo cellum in the low-G+C Gram-positive 
bacteria, recovered from the gut of the termite R. speratus. Clostridium 
termitidis [43] and Termitobacter aceticus were from the termite N. lujae. 
See the legend of Fig. 3 for explanation. 



Fig. 7. Phylogenetic tree showing the positions of 16S rDNA phylotypes 
from the gut of the termite R. speratus, affiliated with clusters comprised 
only of sequences from uncultured Clostridia in the low-G+C Gram-pos- 
itive bacteria. The apparent identity between the phylotypes Rs-J36 and 
Rs-K36 was due to the exclusion of sites containing gaps. See the 
legend of Fig. 3 for explanation. 


were from termite guts (the BCf series, described above), 
and those to 46 phylotypes were from the intestines or 
feces of other animals such as pigs [21], cows [26,27], 
chickens [28] and humans [29], The cluster shown in 
Fig. 6 contained no described bacteria and corresponded 
to the cluster of bovine ruminal clones and Cluster-E of 
Leser et al. [21], 

Three phylotypes, each represented by one to four 
clones, were affiliated with lactic acid bacterial clusters 
(Fig. 9A). Rs-D42 formed a monophyletic cluster with 
Enterococcus, and Rs-C68 and Rs-B70 clustered with Lac- 
tococcus. The closest sequence in databases to Rs-B70 was 
BCf2-19 (94.7% similarity, Fig. 9A) from C. formosanus. 
Seven phylotypes affiliated with Mycoplasma constituted a 
monophyletic cluster with the sequences from C. formosa- 
nus and a sequence amplified front a cell of Koruga bonita, 
a symbiotic protist in the gut of the termite Mastotermes 
darwiniensis [30], as shown in Fig. 9B. The sequence sim- 
ilarities among these termite-derived mycoplasmas were 
79-95%. Three phylotypes represented by one to seven 
clones were affiliated with the Anaerobaculum thermoterr- 
num cluster (Fig. 9C). Two of the three formed a mono- 
phyletic cluster with the two sequences from C. formosa- 
nus, and their similarities were 93-96%. 
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Fig. 8. Phylogenetic tree showing the positions of 16S rDNA phylotypes 
related to C. propionicum (A) and C. purinolyticum (B) in the low-G+C 
Gram-positive bacteria, recovered from the gut of the termite R spera- 
tus. See the legend of Fig. 3 for explanation. 


3.2.4. Termite Group I 

This group of bacteria was first recognized in our pre- 
vious study, as the 16S rDNA clones TG13 (AB004575) 
and TG14 (AB004576) amplified from the gut of R. spe- 
ratus [6], In the present study, five phylotypes were affili- 
ated with this group, and the two phylotypes Rs-D95 and 
R.S-D17 showed high sequence similarities to TGI 3 
(98.3%) and TG14 (99.9%), respectively (Fig. 10). As in 
our previous study [6], clones belonging to this group were 
found abundantly in the sequenced clones, forming one of 


the second most dominant groups (Figs. 1 and 4). The 
similarities between these five phylotypes were more than 
95% and constituted a robust monophyletic cluster (Fig. 

10) . This cluster was related to four environmental sam- 
ples, 4COd-3 from a bovine rumen [27], SBR2095 from 
activated sludge [31], vadinHB76 from wine vinasse [32], 
and a short sequence WCHA1-89 from a contaminated 
aquifer [33], Whereas the high bootstrap confidence value 
(99%) supported the monophyly of the termite cluster and 
the four environmental samples, the sequence similarities 
between them were below 82%. No information, except for 
16S rDNA sequences, has been reported so far on this 
group. 

3.2.5. Proteobacteria 

Thirty-one phylotypes were affiliated with Proteobacte- 
ria and further classified into five subdivisions: 16 a-Pro- 
teobacteria, two y-Proteo bacteria, four (3 -Proteobacteria, 
three e-Proteobacteria, and six 5-Proteobacteria (Fig. 

11) . The 16 a-proteobacterial phylotypes were represented 





Fig. 9. Phylogenetic tree showing the positions of 16S rDNA phylotypes 
affiliated with Bacilli (A), Mycoplasma (B) and Anaerobaculum thermo- 
terrnum cluster (C) in the low-G+C Gram-positive bacteria, recovered 
from the gut of the termite R. speratus. The isolates Enterococcus sp. 
RfL6 and Lactococcus sp. Tml05 were from the guts of the termites 
R. Jlavipes and Thoracotermes macrothorax, respectively [44]. The clone 
Mycoplasma sp. (AJ 132469) is a short sequence (corresponding to posi- 
tions 255-1388 in E. coli JO 1695), amplified from a single cell of the 
protist Koruga bonita in the gut of the termite Mastotermes darwiniensis 
[30]. See the legend of Fig. 3 for explanation. 
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Fig. 10. Phylogenetic tree showing the positions of 16S rDNA phylo- 
types affiliated with the candidate division TG 1 and Acidobacteria 
group, recovered from the gut of the termite R. speratus. The clones 
TGI 3 and TGI 4 were previously recovered from R. speratus [6]. The 
clone WCHA1-89 from a contaminated aquifer [33] is a short sequence, 
corresponding to positions 534-1388 in E. coli (J01695). The phylotype 
Rs-M47 did not cluster with any other sequences in the analysis. See 
the legend of Fig. 3 for explanation. 


each by one to nine clones and the sequence similarities 
were 78-96%, when excluding the phylotype Rs-B62, 
which possessed an about 150-bp insertion in the 5' re- 
gion. All the a-proteobacterial phylotypes were found for 
the first time except for two, Rs-A44 and Rs-K73, which 
showed high sequence similarities to Methylobacterium sp. 
V3 (99.5%) and Methylocella palustris (97.2%), respectively 
(Fig. 11). 

Two phylotypes represented by only one and two clones 
were affiliated with y-Proteobacteria, and the phylotypes 
Rs-N33 and Rs-M74 were almost identical to the sequen- 
ces from Pseudomonas aeruginosa (100%) and Enterobacter 
anmigenus (99.4%), respectively (Fig. 11). Among the four 
(3-proteobacterial phylotypes, Rs-A34 was the most abun- 
dant from the sequenced clones (Fig. 4), and affiliated with 
the Rhodocyclus group (Fig. 11). In the Rhodocyclus 
group, the phylotype Rs-K54 shared 95.6% sequence iden- 
tity with BCf6-08 (AB062831) front C. formosanus. Rs- 
D51 was found to be related to Chromobacterium viola- 
ceum (M22510) and almost identical to the sequence 
OS1L-24 (AB086875) (99.5%) from activated sludge. The 
three e-proteobacterial phylotypes were represented each 
by one to five clones and formed a monophyletic cluster, 
showing sequence similarities of 93-95% (Fig. 11). There 
was no cultured bacterium showing a high sequence sim- 
ilarity to them. The 5-proteobacterial phylotypes were 
comprised of five Desulfovibrio and one related to Geo- 
bacter (Fig. 11). Among the Desulfovibrio phylotypes, 
Rs-N31 was the most abundant (Fig. 4), and is related 
to Desulfovibrio intestinalis isolated from the gut of the 
termite M. darwiniensis [34], although the sequence simi- 
larity was only 88%. 

3.2.6. Actinobacteria 

Fourteen phylotypes, represented by one to three clones, 


were affiliated with Actinobacteria (Fig. 12). Five phylo- 
types were clustered with Atopobium minutum (S44204) 
and BCf9-ll (AB062844) from C. formosanus. The other 
five phylotypes were clustered with RFS10 (AF068463), 
part of a chimeric sequence from the termite R. flavipes 
[9], They exclusively constituted a robust monophyletic 
cluster with sequence similarities of more than 88% (Fig. 
12). Rs-N96 was the closest to RFS10, sharing 95.9% se- 
quence identity. 



Fig. 11. Phylogenetic tree showing the positions of 16S rDNA phylo- 
types affiliated with Proteobacteria, recovered from the gut of the ter- 
mite R. speratus. Desulfovibrio intestinalis [34] and D. termitidis were iso- 
lated from the termite M. darwiniensis. Sinorhizobium sp. M3A was 
isolated from the termite Nasutitermes nigriceps. See the legend of Fig. 3 
for explanation. 
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3.2.7. Phylogenetic groups with few representative clones 

The phylotype Rs-M39 was clustered with Acidobacte- 
rium capsulatum (D26171) and shared 97.0% sequence 
identity with the clone TRB82 (AF047646), recovered 
from an acid mine drainage site (Fig. 10). Rs-M47, repre- 
sented by only one clone, did not cluster with any other 
sequences in the analysis (Fig. 10). The sequence showing 
the highest homology to Rs-M47 in the BLAST analysis 
was sequence BD4-10 (AB015569), recovered from deep 
sea sediment [35], although the similarity was as low as 
79%. 

Four phylotypes, affiliated with Planctomycetes, were 
clustered with BCf2-25 (AB062813) from C. formosanus, 
constituting a robust monophyletic cluster (Fig. 13 A). The 
sequence similarities within this cluster were 89-96%. Two 
phylotypes were affiliated with Verrucomicrobia, and Rs- 
D37 was clustered with BCf3-05 (AB062814) from C. for- 
mosanus, with a high bootstrap confidence value (Fig. 
13B), sharing 82.5% sequence identity. Two phylotypes 
were affiliated with the candidate division OP 11 [36] 
(Fig. 13C). The phylotype Rs-J47 was clustered with the 
oral clone XI 12 (91% similarity), and Rs-M28 was clus- 
tered with clones BD1-5 [35] and VC1.2-cll3 [37], recov- 
ered from deep sea sediment (89% similarity). The phylo- 
type Rs-FI34 was clustered with Gleobacter violaceus in the 
phylum Cyanobacteria (Fig. 13D), and shared 94.7% se- 
quence identity with the bovine ruminal clone 4COd-2 
[27]. The remaining two phylotypes, Rs-A23 and Rs-J96, 
did not cluster with any known sequences with a high 
bootstrap confidence value (Fig. 13D). The sequences 
listed in the results of the BLAST and FASTA analysis 
of these two phylotypes were from Clostridia, including 
Caloramator sp. JW/JH-1 (AF286863), Caloramator cool- 
haasii (AF104215) and Clostridium thermopalmarium 



Fig. 12. Phylogenetic tree showing the positions of 16S rDNA phylo- 
types affiliated with Actinobacteria, recovered from the gut of the ter- 
mite R. speratus. RFS10 is a short sequence (corresponding to positions 
64-1111 in E. coli JO 1695), recovered as part of a chimeric sequence 
from the termite R. fiavipes [9]. See the legend of Fig. 3 for explanation. 
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Fig. 13. Phylogenetic tree showing the positions of 16S rDNA phylo- 
types affiliated with (A) Planctomycetes, (B) Verrucomicrobia, (C) the 
candidate division OP 11 and (D) Cyanobacteria, recovered from the gut 
of the termite R. speratus. The ‘unidentified cluster’ was not related to 
any known sequences with a high bootstrap confidence value. See the 
legend of Fig. 3 for explanation. 


(X72869), while the sequence similarities were lower than 
80%. 

In total, 248 of the 268 phylotypes were found for the 
first time in the present study. Some of these phylotypes 
were clustered with each other and/or with sequences de- 
rived from other termite species. This implies the existence 
of termite-specific bacterial lineages. However, it is cur- 
rently impossible to examine the coevolutionary relation- 
ship between termites and the bacterial symbionts because 
only small amounts of sequence data from other termite 
species are available. We expect that this study becomes 
the basis for further studies on the diversity of termite gut 
microbiota and the phylogenetic relationships of these 
bacteria, as well as the localization and characterization 
of the bacteria. 
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